Not to appear in Nonlearned J., 45. 



Adaptive Optics Assisted Near-Infrared Spectroscopy of SVS 13 

and its Jet^ 

C.J. Davis 

Joint Astronomy Centre, 660 N. A'ohdku Place, Hilo, Hawaii 96720, USA 

c . davis® j ach . hawaii . edu 
B. Nisini 

INAF- Osservatorio Astronomico di Roma, Via di Frascati 33, I-OOO4O Monteporzio 

Catone, Italy 

M. Takami and T.-S. Pyo 

Subaru Telescope, 650 N. A'ohdku Place, Hilo, Hawaii 96720, USA 

M.D. Smith 

Armagh Observatory, College Hill, Armagh, N. Ireland 
E. Whelan and T.P. Ray 
Dublin Inst. Advanced Studies, 5 Merrion Square, Dublin 2, Ireland 

A. Chrysostomou 

Centre for Astrophysics Research, University of Hertfordshire, Hatfield, Herts ALIO 9AB, 

UK 

ABSTRACT 

We present long-slit H- and K-band spectroscopy of the low-mass outflow 
source SVS 13, obtained with the adaptive-optics assisted imager-spectrometer 
NACO on the VLT. With a spatial resolution of <0.25" and a pixel scale of 



"'^Based on Observations collected at the European Southern Observatory, Paranal, Chile (ESO pro- 
grammes 074.C-0408). 
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0.027" we precisely establish the relative offsets of H2, [Fe ll], CO, Hi and Nal 
components from the source continuum. The H2 and [Fe ll] peaks are clearly 
associated with the jet, while the CO, Hi and Nal peaks arc spatially unresolved 
and coincident with the source, as is expected for emission associated with accre- 
tion processes. The H2 profile along the slit is resolved into multiple components, 
which increase in size though decrease in intensity with distance from the source. 
This trend might be consistent with thermal expansion of packets of gas ejected 
during periods of increased accretion activity. Indeed, for the brightest com- 
ponent nearest the source, proper motion measurements indicate a tangential 
velocity of 0.028"/year. It therefore seems unlikely that this emission peak is 
associated with a stationary zone of warm gas at the base of the jet. However, 
the same can not be said for the [Fe II] peak, for which we see no evidence for 
motion downwind, even though radial velocity measurements indicate that the 
emission is associated with higher jet velocities. We postulate that the [Fe II] 
could be associated with a coUimation shock at the base of the jet. 

Subject headings: ISM: jets and outflows - ISM: Herbig-Haro Objects - stars: 
pre-main-sequence - stars: individual (SVS 13) 

1. Introduction 

With high-spectral-resolution observations one may probe emission-line regions at the 
bases of jets, and thereby hope to distinguish between - or at least constrain - models of jet 
coUimation and acceleration. Indeed, there are a number of diagnostic atomic and molecular 
lines in the near-infrared (near-IR) that are particularly useful for studying the more deeply 
embedded Young Stellar Objects (YSOs) (Davis et al. 2001, 2003; Pyo et al. 2002; Whelan 
et al. 2004; Takami ct al. 2005; Nisini ct al. 2005). Most notable are the [Fe ll] and H2 
lines, which are found to be bright at the base of a number of jets from Class 1 protostars. 
The emission is usually blue-shifted, the [Fell] tracing a higher-velocity component than 
the H2, and offset along the axis of the jet by a few 10s of AU. These two species trace 
very different ffow components: the [Fe ll] derives from hot, dense, partially- ionised, gas 
(T~10,000K) while the H2 traces low-excitation, shocked (or possibly fluoresced) molecular 
gas (T~2,000K). These forbidden emission line (FEL) and molecular hydrogen emission fine 
(MHEL) regions exist within a few hundred AU of each outflow source, and are therefore 
either within or very close to the primary jet coUimation and acceleration zone. The H2 may 
be entrained along the walls of the [Fe ll] jet, or it may trace a cooler molecular disk wind 
ejected at larger disk radii. 
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As a further step toward a better understanding of these emission-hne regions, here 
we present adaptive-optics (AO) corrected H- and K-band spectra of the MHEL and FEL 
regions associated with SVS 13, the driving source of the HH 7-11 Herbig-Haro (HH) jet and 
molecular outflow (Chrysostomou et al. 2000; Bachiller et al. 2000). The source is located 
in the L 1450 molecular cloud near the active low-mass star forming region NGC 1333 in 
Perseus. SVS 13 is assumed to be at a distance of 300 pc (Cernis 1990; de Zeeuw et al. 
1999); modeling of the HH bow shocks in the flow imply an inchnation angle of ~30°-50° 
with respect to the plane of the sky (Hartigan, Raymond & Hartmann 1987; Smith, Khan- 
zadyan & Davis 2003). The FLT observations presented here complement the earher U.K. 
Infrared Telescope (UKIRT) observations of Davis et al. (2001, 2003), which yielded impor- 
tant kinematic information though lacked spatial resolution, and the more recent Subaru 
observations of Takami et al. (2005), where a thorough excitation analysis has been made. 



2. Observations 

The combination of NAOS and CONIC A at the European Southern Observatory (ESO) 
Very Large Telescope (VLT) provides AO correction using an IR- bright guide star in conjunc- 
tion with near-IR spectroscopy. In our case, the jet source itself, SVS 13 (J~11.6, K~8.1), 
was used as the AO reference star. The S27-3-SH and S27-3-SK modes used yield a wave- 
length coverage of 1.37-1.72/im and 2.02-2. 53yum at a spectral resolution of R~1500. The 
slit width was 0.172", the pixel scale along the slit was 0.0270"(±0.0002") (measured from 
the offsets of a bright star along the slit), and the slit length was 28". The H- and K-band 
observations were conducted in service mode on 13 and 17 January 2005, respectively. With 
the slit aligned along the SVS 13 jet axis (position angle = 159° E of N), the telescope was 
nodded between object and blank sky six times in each waveband. Two 150 sec exposures 
were obtained at each position. 

The individual spectral images were bad-pixel masked and flat-fielded using observations 
of a halogen lamp. Sky exposures were subtracted from object frames before wavelength 
calibration using argon arc spectra. Starlink FIGARO routines were used to correct for 
distortion along the slit axis, so that the wavelength calibration was constant along the spatial 
axis (along columns). The continuum associated with SVS 13 was then registered along the 
slit axis before co-addition of the data. Individual sky-subtracted images were compared to 
the shifted and co-added data in each waveband to make sure that the spatial resolution had 
not been compromised by this shifting and averaging process. From Lorentzian fitting to 
profiles taken through the SVS 13 continuum at various wavelengths, the full width at half 
maximum (FWHM) was found to vary between 0.24" and 0.40" in H, and 0.19" and 0.35" 
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in K; in the averaged data the FWHM was 0.26"(±0.03") in H and 0.23"(±0.02") in K. 

The B5V bright star HIP 30943 (V=8.0) was observed as a telluric standard. However, 
the extracted H- and K-band spectra of this source contain deep Hi brackett absorption 
hues. Model Vega spectra were therefore used to remove these hues. Briefly, via an iterative 
process similar to that described by Vacca, Gushing & Rayner (2003), the model data were 
scaled, smoothed and shifted before division into the HIP 30943 spectra. Comparison of the 
corrected HIP 30943 spectra with theoretical atmospheric transmission spectra presented by 
Vacca et al. (2003) showed that this process worked very well; the small difference in spectral 
type between Vega (AOV) and the standard (B5V) had no noticeable affect on the Hi fitting. 
The atmospheric absorption features remaining in the corrected HIP 30943 spectra were 
subsequently aligned with equivalent features in the H- and K-band SVS 13 data, before the 
HIP 30943 spectra were "grown" into spectral images and divided into the SVS 13 spectral 
images. 



3. Near-IR spectroscopy 

H and K-band spectral images are shown in Figures 1, 2 and 3. To more clearly show 
the line-emission features, the SVS 13 continuum has been removed from each image, row 
by row, by fitting a second-order polynomial to wavelength regions free of line emission. 
However, residual shot noise associated with the bright continuum does remain in some 
parts of the data. 

As is the case in other spectroscopic studies (Davis et al. 2001; Takami et al. 2005), 
extended emission is observed only in H2. However, for the first time, faint H2 emission is 
also observed in the SVS 13 counter-jet (positive offsets in Figure 3). The HI Brackett fines 
in the H-band, and the CO, Br7 and Nal lines in the K-band all appear as compact peaks 
coincident with the continuum. The [Fe ll] peak is also compact, although it is slightly offset 
along the blue-shifted jet. 

In Figure 4 we compare plots of the H2 l-OS(l) and [Fe ll] 1.644 //m profiles traced 
along the slit, together with profiles of the adjacent continuum. The continuum plots are the 
average of profiles extracted slightly blue-ward and red-ward of the line in each case. The 
continuum profiles effectively show the spatial resolution of the observations, although they 
may be broadened slightly by nebulosity associated with SVS 13. Similar profile plots were 
produced for the CO 2-0 bandhead and Br7 line emission (not shown). 

Only one component is evident in each of the CO, Br7 and [Fe II] profiles, while as many 
as five separate peaks are identified in the H2 profile (labeled in Figure 4). From Lorentzian 
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fitting, the offsets with respect to the source continuum position and the FWHM of the 
individual components seen in each line were measured. These are given in Table 1. The 
CO and Bvy peaks are unresolved and are precisely coincident, to within a few AU, with the 
source continuum; we see no evidence for CO or Hi components associated with the outflow. 

The [Fe ll] profile is, on the other hand, marginally resolved. In the echelle data of Davis 
et al. (2003) and Takami et al. (2005) a radial velocity of ~ 140 km s"-*^ with respect to the 
systemic velocity was assigned to this component. Presumably this compact, high- velocity 
emission derives from a discrete knot or shock front near the base of the jet, rather than 
from an extended emission region. 

In the H2 data, we see an interesting trend where the size (FWHM) of each component 
increases with distance from the source (Table 1). Moreover, if the H2 and [Fe II] trace the 
same jet material, this trend extends to the [Fe II] data; the closer a component is to the 
source continuum, the less extended it appears to be. The H2 components also decrease in 
intensity with distance. This behavior could be explained by the coohng, decreased excitation 
and expansion of "packets" of gas as they travel away from the source. A decrease in gas 
density and excitation temperature is inferred from the excitation analysis of Takami et al. 
(2005). A similar trend is also seen in small-scale jets from T Tauri stars (Bacciotti et al. 
2000; Woitas et al. 2002). 

Directly behind a fast- moving shock the very hot (T ~ 10^ — 10^ K) post shock gas will 
radiatively cool very rapidly, before the region has time to expand. However, the bulk of the 
material in the gas packet, behind the shock front working surface, may expand adiabatically. 
If the pressure inside the warm (T ~ 100 K) gas packet approaches that of the surrounding 
medium as it expands, decreasing by a factor of ~10-100, and if at the same time the 
gas temperature in the packet decreases by a factor of 3-10 as it travels a distance of a 
few arcscconds (equivalent to the inter-knot spacing), then the packet might be expected to 
increase in volume by a factor of 1-30, or in diameter by a factor of 1-3. Isothermal expansion 
would result in greater expansion, while a lower pressure gradient would suppress expansion. 
Either way, the ratio is potentially consistent with the relative sizes of the H2 components 
in Table 1. 



4. Proper Motions 

By comparing the VLT data with the earlier Subaru observations of Takami et al. (2005) 
we can measure, or at least set upper limits on, the proper motions (PMs) of the [Fe ll] and H2 
components. The .S'u&aru observations, although at higher spectral resolution (i? ~ 1.1 x 10^), 
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were not obtained with AO correction. However, the Subaru pixel scale (0.060"(±0.002")) 
does fully sample the seeing and the same slit position angle (159°) was used with both 
instruments. Offsets of emission features along the jet axis can therefore be measured accu- 
rately and a direct comparison made between the VLT and Subaru observations (sec Table 2). 
Note, however, that the Subaru slit was slightly wider than the VLT slit; 0.30" versus 0.172". 

• In the VLT H-band observations, the single [Fe ll] component in Figure 4 may cor- 
respond to the high-velocity peak seen in the Subaru data, for which Takami et al. 
(2005) assign a radial velocity of ~ 140 km s~^. Gaussian fitting yields an offset of 
0.10"(±0.01") for the Subaru peak. Curiously, the offset of the FLT peak in Table 2 
implies that this feature has moved upwind, i.e. closer to the source. We certainly see 
no evidence for pronounced movement away from the jet source. 

• In the Subaru K-band echelle data two velocity components are resolved in H2; a bright, 
low-velocity component (LVC) blue-shifted by ~30 km s~^ that is spatially offset along 
the jet by 0.23"(±0.01"), and a fainter, more diffuse, high- velocity component (HVC) 
blue-shifted by ~100 km s"^ that is offset by 1.27"(±0.02"). These two peaks, which 
are spatially as well as kinematically separate, probably correspond to the first two 
components in Figure 4. In Figure 5 we compare the Subaru data with smoothed VLT 
data: here the components clearly have a similar spatial offset, width and brightness 
ratio. (Note that if the data were normalized to the HVC/component 2 peak, then 
the Subaru LVC would be brighter than the VLT component 1. Given the wider 
Subaru sht, one would expect this to be the case if the LVC traced a broader flow 
component than the HVC.) For the LVC/component 1, the offsets in Table 2 indicate 
a shift downwind of 0.06"(±0.02") over the 783 day time interval between the Subaru 
and VTT observations. However, the HVC/component 2, like the high- velocity [Fe ll] 
component, again exhibits no measurable PM away from the source. 

From the radial velocities of the HVC and LVC, and assuming an inclination angle of 40° 
to the line of sight and a distance of 300 pc to SVS 13, then over the time period between the 
observations the H2 LVC and HVC should have moved by 0.04" (12 AU) and 0.13" (38 AU) 
on the sky, while the [Fe ll] feature should have moved by 0.18" (53 AU). Clearly, only the 
H2 LVC has a PM consistent with the observed radial velocity. Such consistency is expected 
if, for example, the gas is accelerated by a fast-moving shock front in a heavy jet. 

The apparent lack of movement of the H2 HVC could result from: (1) the larger errors 
associated with the positional measurements of this fainter, more extended feature; (2) 
blending of component 2 with component 3 (Table 1) in the lower-spatial-resolution Subaru 
data, so that the overall peak appears shifted downwind; or (3) the fact that the time 
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interval between the observations is comparable to the molecular cooling time and therefore 
the timescale for morphological change, which could introduce unknown errors on the PM 
measurements. 

The lack of a measurable PM in the [Fe ll] peak is more difficult to explain. If the [Fe 
ll] peak was associated with ejected clumps or bullets, then one would perhaps expect to 
see additional components along the flow axis (as is the case in H2). Instead, the [Fe ll] is 
confined to a single, very compact peak at the jet base. The [Fe II] could be associated with 
a stationary, collimating shock, similar to that described in the models of Ouyed & Pudritz 
(1994). If this is indeed the case, the coUimation point would be at a distance of ~ 20 AU 
from the source. 

On the other hand, we may be observing two completely independent [Fe II] peaks in 
the Subaru and VLT data, rather than the same near-stationary feature. The Subaru feature 
may have faded, to be replaced by a new peak in the VLT observations. The cooling time, 
from 20,000 K to 7,000 K for dense (10^ - 10^ cm-^), post-shock gas will be of the order 
of weeks or even days (Smith 2003), so morphological changes are certainly possible. This 
might explain the apparent "upwind" movement of this feature. 

A number of groups have attempted to image the emission-line region at the base of 
the SVS 13 outflow. Davis et al. (2002) used a Fabry-Perot (FP) etalon to "boost" the 
line/continuum ratio in their data. The H2 emission at the jet base is certainly evident 
in their image, although the poor spatial resolution of their observations (~1") somewhat 
limits their ability to extract spatial information. They do collapse their image along an axis 
perpendicular to the jet and attempt to subtract the continuum emission from the resulting 
profile (their Figure 4). They identify a "peak" and a "plateau" along the jet axis which 
may correspond to H2 components 2 and 3 in the VLT data. We list the offsets of these 
features in Table 2. 

Noriega- Crespo et al. (2002) used the HST to obtain high-spatial resolution near-IR 
images of the HH7-11 outfiow. They present a continuum-subtracted II2 image of SVS 13 
(their Figure 3). However, residual artifacts left over from the continuum subtraction, in 
the bright core but also in a ring of radius ~1", make identifying features within 0.5-1.0" 
of SVS 13 more-or-less impossible. However, they do detect and resolve the H2 emission 
in the flow at offsets of >1". Within l"-2" of SVS 13 the H2 appears to comprise at least 
two knots superimposed on to a patch of moderately extended emission, suggesting a broad 
opening angle for the H2 flow. Component 2 in the VLT data (the HVC in the Subaru 
observations) may therefore constitutes two knots; we give a combined offset for these two 
features (measured along the VLT/Subaru slit axis) in Table 2 for reference. 
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Overall, the results from these two imaging surveys can only be used to confirm the 
presence of emission features. Uncertainties in the continuum subtraction limit their useful- 
ness when measuring PMs. Notably, the compact LVC evident in both the Subaru and VLT 
spectroscopy was not extracted from the FP or HST images. 

Clearly, follow-up observations with the NACO system are needed to better constrain the 
PM measurements discussed above, particularly for the [Fe II] peak and the H2 components 
further downwind. Confirmation of the "stationary" [Fe II] peak would add credence to 
the collimation-shock interpretation. If we also ultimately find that the H2 peaks do not 
move, then in H2 we could be witnessing a warm, stationary region or perhaps a turbulent 
mixing layer between the jet and ambient medium at the jet base through which jet material 
flows, rather than discrete bullets or clumps ejected during episodes of increased accretion. 
However, given the tentative result for the LVC, and the fact that knots in Herbig-Haro jets 
and molecular flows typically do have large proper motions (Reipurth & Bally 2001), for H2 
at least, this seems unlikely. 

Finally, we mention that the dynamical age of the H2 LVC, components 1 in the VLT 
data, when derived from its PM and offset in Table 1, is ~ lOyrs. Notably, the optical 
outburst reported by Eisloffel et al. (1991) occurred roughly 15yrs before the observations 
reported here. Our results arc thus reasonably consistent with the H2 LVC being produced 
by this outburst and the associated increase in accretion. 
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Fig. 1. — Contour plot extracted from the K-band spectral image showing the H2 l-OS(l) 
emission. Negative offsets are along the south-eastern, blue-shifted jet axis. The continuum 
associated with SVS 13 has been fitted and removed (see text). 



Table 1: Parameters from the profile fitting 



Line^ 


Offset 


Offset 


FWHM 




(arcsec) 


(AU) 


(arcsec) 


H2 - 1 


0.29(±0.009) 


87.8(±3) 


0.44 


H2 - 2 


1.19(±0.018) 


358(±6) 


1.3 


H2- 3 


2.0(±0.1) 


600(±30) 


1.5 


H2-4 


4.8(±0.2) 


1428(±30) 


1.8 


H2- 5b 


0.28(±0.1) 


86(±30) 


1.0 


[Fell] 


0.061(±0.004) 


18.3(±1.2) 


0.29 


Br7 


0.001(±0.003) 


<0.9 


0.19 


CO 


0.004(±0.003) 


<1.2 


0.17 



"Lorcntzian fits to five components in H2 (see Figure 4); a single Lorentzian is fitted to the 

single peak in each of the extracted CO, [Fe ll] and Br7 profiles. 

^Possible component in the red-shifted counterfiow (positive offsets in Figure 4) . 




Fig. 2. — Contour plots extracted from the H- and K-band spectral images showing emission 
in select lines. The continuum has again been removed. 
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Fig. 3. — Continuum-subtracted K-band spectral image with a grey-scale stretch set to show 
the faint H2 emission along the slit. 
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Fig. 4. — Profile plots showing the distribution of H2 l-OS(l) and [Fe II] 1.644 /xm emission 
along the slit. The fine grey lines show Lorentzian fits to the data; one fit for the [Fe II] peak 
though five fits to the H2 data. The dashed line shows the profile through the continuum 
adjacent to each line. 



Camiurlsan of tha Subani ind VLT H2 prallis 




-ZJ -3 -l.i -1 -DLS D 0.5 1 U 
Qfftatfansci} 



Fig. 5. — Comparison of the Subaru (Grey) and VLT (black) H2 profiles. The dashed lines 
show the adjacent continuum data. The VLT observations have been Gaussian-smoothed to 
the spatial resolution of the Subaru data. 
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Table 2: Component offsets and proper motions 



Line/component^ 




Subaru'^ 


UKIRT^ 


HST 


P.M/ 




Offset (") 


Offset (") 


Offset (") 


Offset (") 


(7yr) 


H2 - 1 


0.29(±0.009) 


0.23(±0.01) 






0.028(±0.007) 


H2 - 2 


1.19(±0.018) 


1.27(±0.02) 


1.28(±0.01) 


-l.l(±0.5)s 


<0.05 


H2 - 3 


2.0(±0.1) 




2.13(±0.09) 




<0.03 


H2-4 


4.8(±0.2) 










[Fell] 


0.061(±0.004) 


0.10(±0.01) 






-0.018(±0.006) 



"Same as Table 1. 

^VLT spectroscopy presented in tliis paper; data obtained on 1/17/2005. 

'^Suharu spectroscopy presented in Takami et al. (2005); data obtained on 11/25/2002. 

'^C/i^/i?r Fabry- Perot imaging from Davis et al. (2002); data obtained on 11/23/2000. 

^HST-NICMOS \m.&g\iig from Noriega-Crespo et al. (2002); data obtained on 1/9/1998. 

■^Proper motion measurements or upper limits; the negative value for the [Fe ll] 1.644 jim. 

peak reflects its apparent motion towards the source. 

^Multiple components spread over an area of roughly 1". 



